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Abstract 

In wave turbulence, it has been believed that statistical properties are well de- 
scribed by the weak turbulence theory, in which nonlinear interactions among 
wavenumbers are assumed to be small. In the weak turbulence theory, sepa- 
ration of linear and nonlinear time scales derived from the weak nonlinearity 
is also assumed. However, the separation of the time scales is often violated 
even in weak turbulent systems where the nonlinear interactions are actually 
weak. To get rid of this inconsistency, closed equations are derived without as- 
suming the separation of the time scales in accordance with Direct-Interaction 
Approximation (DIA), which has been successfully applied to Navier-Stokes 
turbulence. The kinetic equation of the weak turbulence theory is recovered 
from the DIA equations if the weak nonlinearity is assumed as an additional 
assumption. It suggests that the DIA equations is a natural extension of the 
conventional kinetic equation to not-necessarily-weak wave turbulence. 

Keywords: wave turbulence, turbulent statistics, direct-interaction 
approximation 



1. Introduction 

Energy in wave fields such as ocean surface waves [l| , ocean internal waves Q , 
and elastic waves on thin metal plates Q is transferred among wavenumbers ow- 
ing to nonlinear interactions. The large degree-of-freedom wave fields are called 
wave turbulence. The weak turbulence theory, in which the nonlinear interac- 
tions are assumed to be small, is partially successful in wave turbulent statistics. 
Therefore, "the wave turbulence" and "the weak turbulence" are often regarded 
as synonyms. 

However, the separation of the linear and nonlinear time scales, which is 
assumed in the weak turbulence theory, is often violated as pointed out based on 
the weak turbulence theory 0, S 0] ■ Especially in anisotropic wave turbulence 
such as ocean internal waves [2| and Alfven waves Q, the separation of the 
time scales is almost always violated. The violations are often observed also in 
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both observations and direct numerical simulations in ocean surface waves [8| 
and ocean internal waves @. In the direct numerical simulations, the violations 
appear to be caused by fast non-resonant interactions. Therefore, the wave 
turbulence theory applicable to the non-weak wave turbulence is required for 
statistics of the not-necessarily-weak wave turbulence. 

In Navier-Stokes turbulence which has strong nonlinear interactions among 
wavenumbers, direct-interaction approximation (DIA) (Io| is considered to suc- 
cessfully describe turbulent statistics of Navier-Stokes turbulence. In this Let- 
ter, DIA is applied to a three-wave turbulent system, and a wave turbulence 
theory that accepts short-time and strong nonlinear interactions is constructed. 



2. Weak Turbulence Theory 

Wave turbulent systems with nonlinear interactions among three wavenum- 
bers is generally governed by a canonical equation 

■ Mp) _ sn 

1 &t 5a* (p) { ' 

for a complex amplitude a(p) with Hamiltonian 

H = {2n)- d J dp uj\a{p)\ 2 

+ (2ir)- 2d J dp d Pl dp 2 (T p ° lP2 a*(p)a( Pl )a(p 2 ) + c.c.) , 

where p is a d-dimensional wavenumber vector, and is a frequency of 

wavenumber p given by a linear dispersion relation. A matrix element T p p lP2 
gives strength of three- wave nonlinear interactions among p, pi and p 2 . A 
functional derivative with respect to a is expressed by S/Sa. Moreover, a* 
denotes the complex conjugate of a, and c.c. also denotes the complex conjugate 
of the previous term. A d-dimensional periodic wave field in a volume (27r) d is 
considered for simplicity. Since the wave media are assumed to have spacial 
inversion symmetry, u(p) = w(-p) and 7p P lP2 = TI pi _ P2 - 
The canonical equation ([TJ is rewritten as 

da{p) . , . , . 
—gf- = -iw(p)o(p) 

-iE 7 ? lP2 a (fO«(P2)-2i J2 r p p ^a*( Pl )a(p 2 ). (2) 

P=Pl+P2 P=~Pl+P2 

If the nonlinear terms in Eq. @ is neglected, it can easily be integrated as 

a(p) = A(p) exp(— iw{p)t + iO(p)). 

Therefore, a(p) is called complex amplitude and represents the behaviour of a 
mode p in the phase space. 
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The turbulent statistics governed by Eq. ^ is conventionally given according 



to the weak turbulence theory [ljj, ■ Equation ([2]) is multiplied by a*(p 4 ) 
and added complex conjugate with p and p 4 interchanged, and the ensemble 
averaging (~) is taken: 



dn(p) 



S pp* = ~i Y] T^ lP2 a*{p i )a(p 1 )a(p 2 ) 



q "PP4 'piP2 l 
P=Pl+P2 



-2i Yl ^ 2 1 *«*(P4)a*(pi)a(p 2 ) 

P=-Pl+P2 

+ c.c.(p o p 4 ), (3) 

where c.c.(p -f-> p 4 ) denotes taking complex conjugates with p and p 4 inter- 
changed. Kronecker's delta is denoted as 8 PiPj . The wave action n(p), which is 
energy of wavenumber p divided by ui, is defined by a{pi)a*(pj) — n(pi)S PiPj . 

By applying the random phase approximation in the weak turbulence theory 
in the most primitive sense a(pi)a* (pj)a* (pk) = to the third-order correlation, 
the nonlinear term vanishes and 

dn{p) = 
dt 

Then, the non-zero third-order correlation is obtained by considering the time 
evolution of the third-order correlation. The time evolution of the third-order 
correlation is written as 

\P=P3+P4 P+P3=P4 / 

+ c.c.(p o pi) + c.c.(p o p 2 ) 
= -2\T^ 1P2 {nin 2 - n(n\ + n 2 )) . (4) 
The random phase approximation is applied to the fourth-order correlation: 



a pi a pk a Pj a pi a pi a pi a pj a pk 



= n{pi)n{pj) (S PiPk S PjPl + 6 PtPl 6 Pj 



Pk, 



The frequency difference among the three wavenumbers is denoted by AujP ipo = 
uj(p) - w(pi) - uj(p 2 ). 

The time variation of wave action n appearing in the right-hand side of 
Eq. (jll) is possibly negligible by comparing with 1/AujP iP2 if the nonlinearity 
is weak. Here, the linear timescales are assumed to be much faster than the 
nonlinear time scales. The separation of the time scales is nontrivial and is 
often violated. 

When the separation of the time scales are valid, Eq. Q can be integrated 
from to to to + r, under the initial condition aa^a^to) = 0. The third-order 
correlation is obtained as 

— — _ - 2i 7p P 1 p 2 (™i™2 - n(m + n 2 )) (exp(-iAwP p2 r) - l) 



1 2 — • A P 

-iAw PlP2 



3 



By employing 



i(exp(-iAwr) - 1) / 1 , 

; = P.V. — + i7T(5(Aw) as r -> oo, 



— iAo; \Au} 
Eq. (|4]) finally results in the kinetic equation 



^ = 4ir( \Vi P2 \ 2 K n 2 - n(m + n 2 )) <5(Aw£ ip J 

\p=Pl+P2 

- E |7 p p 2P | 2 (« 2 n-n 1 (n 2 + n)) ( 5(A^ p ) 



Pl=P2+P 



- E l^mnm-^Cn + m))^^) . (5) 

P2=P+Pl / 

The kinetic equation indicates that the energy is transferred among wavenum- 
bers which satisfy the resonant conditions: 

ip = Pl+P2 ipi=P2+P ip2=P + Pl 

1 OJ = OJi + 0J2 ' } U>1 = + ' 1 0J2 = U> + OJi 

3. Direct-Interaction Approximation for Wave Turbulence 

As written in SjH the weak turbulence theory assumes, derived from the 
weak nonlinearity, that the time scales of the nonlinear energy transfer is much 
larger than the linear time scales. However, in almost all the anisotropic wave 
turbulence, the kinetic equation ([5]) itself gives short-time strong nonlinear in- 
teractions, and the separation of the linear and nonlinear time scales are often 
violated. 

In this section, the direct-interaction approximation (DIA) [l3| , in which the 
separation of the time scales is not assumed but the largeness of the degrees of 
freedom is assumed, is applied to the wave turbulence. Note that the largeness 
of the degrees of freedom is also assumed implicitly in the weak turbulence 
theory. Instead of complex amplitude a, variables b\ and 62 defined by 

h(p)= a(p)+ ;* { - p \ Mp)= a(p) ~ 2 f ( ~ p) , 

are employed. For example, the variables are defined as 



where rj(p) and 4>(p) are the Fourier transform of the surface elevation and that 
of the velocity potential in deep capillary waves , and they are 

y/uN y/g\k\ 
Mp) = « n(P) ' &2(P) ^^ (P) ' 
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where II (p) and 4>(p) are the Fourier transform of the stratification thickness 
and that of the velocity potential in ocean internal waves l4 |. Since &i(p) is 
the Fourier transform of real functions, b{(p) = b*(—p). 
The governing equation is rewritten as 

l dt = £ij(p)bj(p) +^2K k{p,Pi,P2)b ] {-pi)b k {-p 2 ). (6) 

P+Pl+P2=0 

The Einstein summation convention is employed. The linear and nonlinear 
matrix elements are 



Zijip) = U - i)^(p), 
~~2 



M ife (p,pi,p 2 ) = ^ ((-i) i+i 7^ pl _ P2 + (-i) j n-p 2 + (-^p-t 



+ (c.C.(p,pi,p 2 ) -> (~P, -Pl,-P2)) , 

where c.c.(p,pi,p 2 ) — > (— p, — Pi, — P2) denotes taking complex conjugates and 
replacing (p, pi, p 2 ) by (— p, — pi, — P2) simultaneously. 

When a perturbation is added to the jth component of a wavenumbcr p' 
at a time t' ', which is bj(p',t'), the ith component of another wavenumber p 
at a later time t, which is bi(p,t), responds to the perturbation. The response 
function is defined as 

&ij{p,t\p,t ) - 



Sbj(p',t') 

From Eq. (j6)), the governing equation of the response function is given by 
dG m (p,t\p',t') 



C ij (p)G jn (p,t\p',t') 

2 ^2 ■^'ijk(P,Pl,P2)bj{-pi)Gkn{~P2,t\p',t'). 



dt 



P+Pl+P2=0 

The initial condition of the response function is given as 

G ij (p,t'\p',t') = SiAp-p')- 

A perturbation that removes a triad interaction among po, <7o and »"o that 
satisfies po + Qo + *"o = is added to a wave field at a time t' . When the 
degrees of freedom is large enough, the effect of the perturbation is small. The 
variable bi is resolved into no direct-interaction field (NDI) where the direct 
interaction among po, <7o and J"o is removed and direct- interaction field (DI). 
Then, bi — of' + b^' and |&j ^| » \b\ \. Similarly, the response function is also 
resolved into G§ } in NDI and G§ } in DI: G tJ = G^ +G$ and |G^ 0) | > \G$\. 

The equation of b^ is given as 

d -^M =Cij {p)bf\p)+ y: K jk (p,pi,p 2 )bf(- P1 )b^(-p 2 ), 

P+Pl+P2=0 

{p,Pi,p 2 }#{po,qo,i-o} 
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and that of is 



+ 2 Y, X ji (p,Pi,Rl)i>f ) (-Pi)it"(-P2) 

P+Pl+P2=0 

{p,Pi,P2}#{po,9oTo} 

+ M, fe (Po,go,r ) (<5 PPo &f )(_q b )6W(_r d ) + S PPo bf ] (q )b { h 0) (r )) 
+ {Po,go,r }. 

Here, {po,Qo, r o} denotes cyclic permutations from (p , 9o ; r o)- The governing 
equation for G^ ) and G^ is also obtained. The solution of is analytically 
obtained 



bW (p) = [ dt' £ (p, t\p', Wnjk(Po,qo, ro) 



x (Vpo^f (-«),062 ,) (-ro,f) + Wf (go^'^Vo,*')) 

+ {Po,go,r } (7) 

under the initial condition b^p (t = t') = 0. The solution of is also obtained 
analytically and given by b<® and G<-°\ 

To investigate wave turbulent statistics, the correlation between 6's defined 

as 



V ij (p,t,t') = b i (p,t)b* j {p,t') 

is considered. Because of Eq. © the simultaneous correlation is governed by 
the following equation: 



dV%J M ~ ~~ = ^im(p)V m jip,t,t) Nimn(p,q,r)bj(-p)b m (-q)b n (-r) 

p+q+r=0 

+ C.C.(iHj). (8) 

The summation of m and n is taken over 1 and 2. DIA that removes the direct 
interactions among p + g + r = 0is applied to the wave field. First, Eq. ([HJ is 
rewritten as 



dV rj {p,t,t) 



dt 



Cim{p)V m j(p,t,t) 



+ £A/- mm (p,9,r) [bf\-p)b^(-q)b^(-r) + bf\-p)b^(-q)b^(-r) 

p+q+r=0 



+&f (- P )b$(-q)bW (-r) + bf\-p)bt ] (-q)bP (_ r ; 
+ c.c.(i*+j), (9) 
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by resolving b in the right-hand side of Eq. (J5) into (NDI) and (DI). In 
NDI field which has no interactions among p, q and r, &^ (— p)6m' ) (— q)bn^ (— r) = 
0, since 6^(p), &m^(?), fri ^'") are statistically independent. Second, the solu- 
tions ([7]) are substituted into b^ {—p)bm\—q)bn (—f) and similar terms. At 
last, Eq. dSJ is expressed by the different-time correlation and the response 
function: 



dVa(p,t,t) 



Cim(p)Vmj(p,t,t) 



2^JV,„(p,q,r) / dt' (jV 6 c (p,qf,r)G? i0 (-p,t|-p,t') T4 m (g,i',f) V cn (r,t',t) 

,+q+r=0 "'*0 



P+<3+ 

+ AA abc ( 9 ,r,p)G mQ (-q,i|-q,i') H„(r,i',i) V C3 {p,t' ,t) 
+Af abc (r,p, q)G na (-r,t\-r,t') V bj (p,t',t) V cm (q,t',t) ' 
+ c.c.(i o j). (10) 

The statistical independence between bi and G mn is assumed. 

A similar procedure can be applied to the different-time correlation. Hence, 
the governing equation of the different-time correlation is obtained as 



dV l3 (p,t,t') 



dt 



= C im (p)V m j(p,t,t') 



+ 2j2^imn(p,q,r) I f dt"M abc (p,q,r)G ja (-p,t'\-p,t") V bm (q,t",t) V cn {r,t",t) 

p+q+r=0 V *° 

+ f dt" (Ar ab c(q,r,p)G ma (-q,t\-q,t") V bn (r,t",t) V cj (p,t",t') 
J to V 

+Mabc(.r,p,q)G na (-r,t\-r,t") V bj (p,t",t') V cm {q,t" ',*))) ■ ( n ) 
Moreover, the governing equation of the response function is also obtained as 



dG in {p,t\p,t>) 



dt 



= £ij(p)GjntP,t\p,t') 



^J2^ijk{P,q,rWabc{r,p,q) f dt"G ka (-r,t\-r,t") G bn {p,t"\p,t>) V cj (q,t",t). 

_ i _ I n Jt' 

(12) 



p+q+r=0 



Equations (fTUHT2l give a closed equation system and they are the DIA equations 
of the wave turbulence. 

It should be emphasized that the weak nonlinearity is not assumed in the 
procedure. Therefore, the DIA equations of the wave turbulence (|101fl"2"|) can be 
applied also to strongly nonlinear wave turbulent systems. 
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4. DIA Equations for Autocorrelation of Complex Amplitude 



In the previous section, DIA is applied to the variable b. Turbulent statistics 
is described with the complex amplitude a to compare the weak turbulence the- 
ory in this section. The correlation of a is defined as a(p, t)a* (p, t') = N(p, t, t') 
and a(p,t)a(—p,t') — M(p,t,t'). The correlation of b, that is Vij is expressed 
by M and N as 



(p,t,t>) = i (i-^N(p,t,t')+i l -iN*(-p,t,t') 

- (r( <+ ^M(p,t J t / )+i <+J M*(-p,t J t / ))) • 

The initial condition of the cross-correlation is that M(p, to, to) = since a(p) 
and a{—p) are uncorrelated initially. The cross-correlation at later time is much 
smaller than the auto-correlation, \M(p,t,t')\ <C \N(p,t,t')\. Therefore, 

%ip^\pJ) = \ (r^- j) N(p,t,t') +i i -*N*(-p,t,t')) , 
Gij(p,t\p,t') = ± (i-^G(p,t,t')+i l ->G*(~p,tX, 
Equation (fTU|) is rewritten to a equation for N(p,t,i) as 
dN(p,t,t) 



2^ fdt'(G*(p,t,t') (\T p q 



J - t , ; , q _ r \ 2 N*(-q,t',t)N*(-r,t',t) 

p+q^ 



+ \T- p i\ 2 N*{-q, t', t)N(r, t', t) + \T pq r \ 2 N(q, t', t)N*(-r, t', t)) 

- N(p, t', t) (\T^ q _ r \ 2 G(-q, t, t')N*(-r, t', t) 

+ \T rp q \ 2 G(-q 7 t,t')N(r,t\t) - \T p - r \ 2 G*(q,t,t')N*(-r,t',t)) 

+{g o r}) +c.c. (13) 

Similarly, the different-time correlation N(p, t, t') is also obtained from Eq. (fTTj) . 
Moreover, Eq. (fT2"j) can be expressed by N and G. Namely, the DIA equa- 
tions (fTUl - [T2")) in the wave turbulence as equations for the correlation Vij of b 
can be rewritten as equations for the correlation N for a. 

Furthermore, the fluctuation-dissipation relation, N(p, t, t') — n(p, t')G(p, t, t'), 
is employed. The simultaneous correlation is written as N(p, t, t) = n(p, t) from 
now on. Hence, the DIA equations (jlOl [T2|) can be rewritten as 



dn(p,t) 



dt 

p+q+ 



2 E / dt ' (\ TP i-r I' *>(-»•> - n(p, t') (n(-q, t') + n(-r, t'))) 

fq+r=0 •''o 



G* (p, t, t')G(-q, t, t')G(-r, t, t') 

- \Tr P q \ 2 (n(p, t')n(r, if) - n(-q, t') (n(p, t') + n(r, t'))) G* (p, t, t')G(-q, t, t')G* (r, t, t') 

-\r pq r \ 2 (n(p,t>)n(q,t>) - n(-r,f) (n(p, t') + n(q, t'))) G*(p,t,t')G*(q,t,t')G(-r,t,t')) 

+ c.c. (14) 
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and 

dG(p,t,t') 



-iu(p)G(p,t,1?) 



dt 

~ 2 E /* dt " G (P> *"> (l^-r-l 2 t") + n(-r, t")) G(-q, t, t")G(-r, t, t") 

- \r r - q \ 2 (n(-q, t") - n(r, t")) G(-q, t, t")G*(r, t, t") 

~\T pq r \ 2 (n(-r, t") - n(q, t")) G*(q, t, t")G(-r, t, t")) . (15) 

Here, to = t' thanks to the causality. Equations (|14[) and fp~5]) are another set 
ol the DIA equations. 

Equation (|14j) in the DIA equations and the conventional kinetic equation 
look quite similar. Nonetheless, there exist some differences. Hysteresis is 
included in the DIA equations since the time integration is incorporated in 
Eq. (fT4|) . It is in contrast with the Markovian properties of the kinetic equa- 
tion. Furthermore, the response functions appear in the integrand of Eq. (fl4|) . 
Hence, the DIA equations accept nonlinear changes of phases. The coherent 
structures such as freak waves which develop abruptly have the non-Markovian 
properties and the phase entrainment. Therefore, it is expected that the inter- 
mittency is evaluated statistically by the DIA equations. 

The quadratic energy J2 P w (p) n (p) is n °t strictly conserved for Eq. (fT4|) since 
Hamiltonian which is the sum of quadratic and cubic energies is conserved. The 
quadratic energy is conserved for the kinetic equation <j5j> owing to the resonant 
interactions. Therefore, this non-conservation laws implies that the fast non- 
resonant interactions associated with strongly nonlinear coherent structures can 
be statistically estimated by the nonlinear parts of the response functions. The 
numerical simulations of Eqs. (fT4|) and (fl5|) are required to evaluate the fast 
non-resonant interactions for the specific wave turbulent systems. The quadratic 
momentum ^2 p pn(p) is a conserved quantity because of the conditions of the 
wavenumbers. Therefore, one can find an equilibrium solution n(p) cx (p ■ 
U) , where U is an arbitrary constant vector. The entropy X)pl°§ n (p) never 
decrease in Eq. (|14[) . In this manner, some of statistical natures are the same 
as the conventional kinetic equation (JS|). 

For the short-time limit t —>■ to, Eq. (|T4"| is rewritten as 

„4(t _ ^ ( J2\T q P r\ 2 (nfa)n(r) - n(p)(n(q) + n(r))) 

\p—q-\-r 

-^2\V P \ 2 (n(r)n{p) - n(q)(n(r) + n{p))) 



q=r+p 



VMtt I 2 



I pi 

r=p+q 



r; q Hp)n(q) - n(r)(n(p) + n(q))) . (16) 



This is also consistent with the short-time kinetic equation in Ref. [12]. Equa- 
tion (|16[) is valid even for strongly nonlinear regimes since the equation is derived 
without employing the separation of the linear and nonlinear time scales. 
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5. Recovery of Kinetic Equation from DIA Equations 



To recover the kinetic equation from Eq. (|13[) . the weak nonlinearity is as- 
sumed as an "extra" assumption in this section. The time variations of the 
different-time correlation and the response function are as follows: 

dN(p,t,t') . Ar/ 

= —iLji\(p,t,t ) + nonlinear terms, 



at 

dG(p,t,t') 
dt 



= —iuiG(p, t, t') + nonlinear terms. 



The nonlinear terms of the different-time correlation and the response function 
do not contribute to the simultaneous correlation at the leading order. Then, 
by neglecting the nonlinear terms, the leading terms of the different-time corre- 
lation and the response function are obtained as 

N(p, t,t')= n(p, f) e - iw (*- t, > , G(p, t, t') = e^*"*') , 

under the assumption of the weak nonlinearity. Therefore, since n(p, t) varies 
much slower than 1/Aw, n(p,t) — n(p,t') — n(p). This n(p) is the very wave 
action in the weak turbulence theory. By substituting the leading terms of the 
different-time correlation and the response function into Eq. ()13[) . we obtain 



-■if^'i ^|7 q p r | 2 (n( 9 )n(r)-n(p)(n(q)+n(r))) e [ 

-'to \p=q+r 



dn(p) 
dt 

- \ T rp\ 2 (n(r)n(p) - n{q){n(r) + n(p))) e-^V'-O 



-J2\ T p q \ 2 (n(p)"(9) - n(r)(n(p) + n(q))) e" 1A < 

r=p+q / 

+ c.c. (17) 
Since the separation of the time scales is assumed, to can be set to — oo. Then, 

J* dt'e iA ^(*-*') = i(P.V. (^L-^ m6(AL%,)). (18) 

Finally, Eq. (|T7|) results in Eq. (J5J) . Namely, the kinetic equation in the weak 
turbulence theory can be recovered by an additional assumption, that is, the 
weak nonlinearity to the DIA equations. 



6. Concluding Remark 

6.1. Discussion 

In the short-time limit, Eq. (| 16[) is consistent with what is derived along 
the weak turbulence theory. The conventional kinetic equation is also recovered 
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from the DIA equations in the weakly nonlinear limit. The nonlinear parts of 
the response function are irrelevant to the time variation of the wave action in 
both limits. 

As pointed out in direct numerical simulations of four-wave weak turbulent 



system 15], the convergence to the kinetic equation is much faster than the 
convergence of the integral to the S function in Eq. (|18l) . During the intermediate 
time, which is longer than the short time and shorter than the convergence, 
the nonlinear parts of the response function will play a role. The role will be 
evaluated for specific systems with numerical simulations. 

6.2. Conclusion 

The closed equation system is developed for the not-necessarily-weak wave 
turbulence statistics according to direct- interaction approximation (DIA). In 
the procedure, the three assumptions below are made: 

• Quantities in the field without the direct interactions are much larger than 
that in the perturbed field under the largeness of the degrees of freedom 
in the wave field. 

• (p) , b^} (q) and b^ (r) are statistically independent in the field without 
the direct interactions among p, q and r. 

• bi and G are statistically independent. 

By developing the equation system without assuming the weak nonlinearity, 
the DIA equations can consistently make statistical description for the wave 
turbulent system that has even short-time and strong nonlinear interactions. 
It provides us an appropriate tool to evaluate the intermittency in the wave 
turbulent system. 

The kinetic equation in the weak turbulence theory is also recovered from 
the DIA equations. This indicates that the framework of DIA which harness 
the largeness of the degrees of freedom is the natural extension of the weak 
turbulence theory. 
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